Anharmonic effects in the A15 compounds induced by sublattice distortions 
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We demonstrate that elastic anomalies and lattice instabilities in the the A15 compounds are 
describable in terms of first-principles LDA electronic structure calculations. We show that at T = 
VaSi, VaGe, and NbaSn are intrinsically unstable against shears with elastic moduli Cn — C12 
and C44, and that the zone center phonons, r2 and ri2, are either unstable or extremely soft. 
We demonstrate that sublattice relaxation (internal strain) effects are key to understanding the 
behavior of the A15 materials. 

PACS numbers: 62.20Dc, 74.25.Jb, 74.25.Kc, and 74.70Ad 
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The A15 materials have been of interest for more than 
40 years due to the relatively high superconducting tran- 
sition temperatures for many of them and the interplay 
between their superconductivity, their structural phase 
transitions, and their temperature-dependent anomalies 
in properties such as elastic constants, phonon spectra, 
electrical resistivity. Knight shift and magnetic suscepti- 
bility. The paradigm interesting A15 materials are 
VsSi and NbsSn, both of which have relatively high 
values and also undergo cubic-to-tetragonal structural 
phase transitions just above (with f > 1 and ^ < 1, 
respectively), and have many other anomalous proper- 
ties as well. Many theoretical models, discussed in Refs. 
[1-3], attempt to explain the unusual properties of the 
A15 materials in terms of models with a sharply peaked 
structure in the electronic density of states (DOS) near 
the Fermi level, Ey- 

Although the past ten years has seen a major focus 
on developing an understanding of the ceramic oxide su- 
perconductors, Q there remain many unanswered fun- 
damental questions related to the unusual properties of 
the A15 materials, including the question of the relation- 
ship between their high and their structural anoma- 
lies/instabilities, and whether or not they are driven di- 
rectly by the electronic distribution using a conventional 
Fermi liquid description of the electronic structure, or 
whether "exotic" excitations are needed to explain the 
properties. Some of these issues may have relevance to 
the high oxides as welL where structural instabilities 
are also known to exist, jj] 

Here we show that using state-of-the-art first-principles 
electronic structure methods to study three of the most 
interesting A15 materials: VaSi, NbaSn, and VsGe. We 
have found several provocative results, including the facts 
that: (1) for all three materials the calculated elastic con- 
stants agree very well with the observed room tempera- 
ture values if sublattice relaxation effects are frozen out; 
(2) when sub-lattice relaxation effects are included, the 
crystals are unstable in the cubic A15 structure with re- 
spect to a tetragonal distortion which stabilizes both the 
elastic constants and the phonon anomalies; (3) the struc- 
tural instability is weakest (in terms of the energetics) in 



VsGe and may not be observable due to effects of crystal 
imperfections; (4) anomalies in several of the zone-center 
phonon branches occur along with the elastic anomalies, 
the branches corresponding to modes where atoms in the 
chains vibrate in the chain direction. 

The electronic structure results reported in this Let- 
ter were generated using the all-electron, full-potential 
linearized-augmented-plane-wave (LAPW) method in 
the local density approximation (LDA). The va- 
lence/conduction band states were treated semi- 
relativistically (no spin-orbit effects) while the core states 
were treated fully relativistically. Care was taken to en- 
sure appropriate k-point convergence for the accuracy 
needed in these studies. 

Usually it is found that the elastic constants, when 
evaluated at the experimental volume, are within 5-10% 
of the experimental values. |^] For a cubic system (such 
as the A15) to be mechanically stable, the bulk modulus 
[B = Cjjj^Cjj. ) and shear elastic constants (Cn — C12 
and C44) must be positive. The overall phonon spectrum 
must be stable as well (positive square frequencies). The 
elastic constants can be extracted using the procedures 
discussed by Mehl et al: ||] (i) calculate the total energy 
as a function of the volume, and then fit the results to 
an equation-of-state to extract the bulk modulus B; 
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FIG. 1. Schematic diagram of ideal (c = a) and tetrag- 
onally distorted (c 7^ a) A15 crystal structures (^6-62). B 
atoms form a body-centered cube, while A atoms form lin- 
ear chains on each cube face. Arrows denote allowed internal 
degrees of freedom (direction of sublattice relaxations) in the 
tetragonally distorted A15 structure. 
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(ii) apply an appropriate strain to the lattice to distort 
the primitive lattice vectors, calculate the total energy 
versus strain S at volume V, and determine the second- 
order elastic constants from the coefficient second order 
in the strain. For example, using an orthorhombic dis- 
tortion, the distortion energy is, 

AE{±S) ^ViCii~Ci2)S^ + 0(6^) . (1) 

We can thus extract Cn — C\i from Eq. (|l|). Similarly, 
we used a monoclinic distortion to extract C44. 

Our calculated total energies of VaSi, VsGe, and 
NbsSn in the A15 structure as a function of volume, were 
used to extract the equilibrium volumes and bulk moduli 
of these materials. The calculated equilibrium volumes 
are underestimated by 4-6%, typical of LDA calculations, 
leading to the calculated equilibrium bulk moduli being 
overestimated compared to experiment by more than ~ 
25%. On the other hand, the calculated bulk moduli at 
the experimental volumes are very close (~ 5%) to those 
of experiment, as shown in Table |. 

Next, we examined the stability of the A15 materials 
under the application of shears corresponding to the elas- 
tic constants Cn — C12 and C44, first without allowing 
the cell-internal sublattice displacements that are allowed 
by the broken symmetries produced by the elastic distor- 
tions. Figure H depicts our calculated distortion energies 
(open squares) as a function of 6^ using orthorhombically 
distorted cells. For a perfectly harmonic crystal, the data 
points would fall on a straight line, with the slopes given 
by Cii — C12. Indeed the calculated points fall quite close 
to straight lines, indicating only small contributions from 
higher order terms [0(5"*)]. The fitted Cn — C12 and C44 
values, shown in Table |, are very close to the experi- 
mental room temperature values, rather than to the low 
temperature values as might be expected for a typical 
material. Since first-principle calculations are supposed 
to correspond to T = 0, not to high temperatures, why 
don't we reproduce the soft elastic constants seen at low 
temperature? 

This puzzle was solved by realizing that as one distorts 
a crystal from a higher to a lower symmetry space group, 
new internal degrees of freedom could be introduced with 
significant physical and computational importance, espe- 
cially for systems with lattice instabilities. We note that 
it is tempting to overlook these sublattice distortions be- 
cause determining their effect introduces a large increase 
in the computational burden, since for each internal dis- 
tortion a set of self-consistent calculations needs to be 
performed to find the equilibrium sublattice distortions 
corresponding to the volume and strain under investiga- 
tion. 

For the orthorhombically distorted A15 cell, the sub- 
lattice displacements correspond to pairs of A atoms 
moving toward or away from each other (rather than fix- 
ing on their symmetry mandated A15 positions), but 
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FIG. 2. Distortion energies for orthorhombically distorted 
A15 structure. Calculations were performed at the experi- 
mental A15 volumes. The open squares denote calculation 
performed with A15-like internal coordinates, while filled cir- 
cles denote calculations with internal coordinates fully re- 
laxed. 

different pairs move independently. Thus, we need to 
add (iii) sublattice relaxation into our above procedures 
for calculating elastic constants. After fully relaxing the 
sublattices, we have found that the relaxed distortion 
energies for the A15 compounds (denoted by filled cir- 
cles in Figure |^) are lowered by a substantial amount, 
which increases as the distortion gets larger. The re- 
laxed distortion energies show a minimum around d ^ 
1-2% and even dip below the energies of the "ideal" cu- 
bic A15 structure. Hence, these three A15 compounds 
arc intrinsically unstable in the cubic structure at T = 0. 
Note that the Cn — C12 instability appears to be weak- 
est in VsGe. Again, the weakest C44 instability occurs in 
VsGe, while the strongest instability occurs for NbaSn, 
as seen in experiment (Table |). 

An alternative way to obtain Cn — C12 is to distort the 
cell tetragonally. This tetragonally distorted cell turns 
out to have the same structure as the observed tetrag- 
onal phase Q for some A15 compounds (Fig. |l|). The 
distortion energies are depicted as a function of distor- 
tion ((5 = f - 1) in Fig. |. If we freeze the atoms in their 
A15-like coordinates, the calculated sublattice-unrelaxed 
distortion energies (denoted by open squares) fall quite 
closely on a parabola. Note that the distortion energy 



2 



contains odd terms of distortion 5 = - — 1 rather than 

a 

being an even function of S as in the case of orthorhom- 
bic and monochnic distortions, so that we needed to per- 
form more first-principles calculations (for ±5) for this 
case. Indeed, we see a slight asymmetry about the <5 = 
axis, indicating a small cubic term [0((5'')]. Least-square 
fitting these data to the form aS'^ + bS^ + cS*, we can ex- 
tract Cii — Ci2 = 3a/V. Such fitted Cn — C12 values are 
within 6% of those shown in Table | using orthorhombi- 
cally distorted cells. 

Next we allowed the internal coordinates to relax in 
a manner consistent with the tetragonal structure, and 
found that the distortion energies were lowered by a large 
amount, showing a double well structure with respect to 
-. In the case of V^Si, the minimum occurs for - > 1, 
while for NbsSn the minimum occurs for - < 1, in 
agreement with experimental observations 1^. For the 
case of VsGe, we find an extremely shallow well and a 
much weaker instability (bordering on the accuracy of 
our LAPW calculations). Therefore, a small amount of 
sample imperfection, such as impurities, vacancies, or 
disorder, or thermal effects, could inhibit the transition 
from occurring in VsGe. Experimentally, no cubic-to- 
tetragonal phase transition has been observed for VsGe, 
although there is some evidence for anomalies in its elas- 
tic constants 0. 

For VsSi and NbaSn, the details of the structural phase 
transitions are very sensitive functions of sample condi- 
tions, with, for instance, the observed ^ ratios, vary- 
ing from 0.9964 to 0.9938 for NbsSn and approximately 
1.0024 for VgSi, 0, while the calculated minimum occurs 
at f about 1.02 and 0.985 for VaSi and NbgSn, respec- 
tively, significantly larger than experiment. Since the 
relaxed distortion energy curves are very small and the 
curves are very flat, a precise determination of the ^ ratio 
is difficult, although most of the quantitative disagree- 
ment between theory and experiment can be ascribed to 
LDA errors that do not affect the conclusions regarding 
the instabilities that we are focusing on. 

We see from these results that the key to understanding 
the elastic constant softening, and most probably other 
anomalous properties in the A15 materials, is the relax- 
ation of the internal degrees of freedom (sublattice re- 
laxation) when there are cell distortions, which leads to 
highly anharmonic behavior in the Al5's, general ideas 
put forward by Testardi. Q At high temperatures, ther- 
mal motion washes away the well structure and atoms 
spend the majority of their time away from the flat 
well minima and assume "average" A15 positions. Cal- 
culations with "frozen-in" Al5-like coordinates give the 
maximum attainable elastic constants, which correspond 
closely to the experimental room temperature values. As 
the temperature is lowered, the effect of sublattice re- 
laxation sets in, and at T = 0, these A15 compounds 
are intrinsically unstable with respect to Cn — C12 and 
C44 shear distortions. These results can be interpreted 



in terms of a ferroelastic type of instability as discussed 
by Miller and Axe ||l^ , with examples given in the book 
of Silje. 

We have seen that the first-principles electronic struc- 
ture method predicts the existence of elastic constant in- 
stabilities and the cubic to tetragonal phase transition in 
A15 materials. It is also well established that the overall 
phonon spectra of many of the A15's exhibits unusual 
temperature dependencies and selected softening in dif- 
ferent parts of the Brillouin Zone |]6|,p|,p^ . Such insta- 
bilities for the phonons also follow from our calculations 
using the frozen phonon method which we have applied to 
determine the zone-center optic phonons. Note that the 
relaxed distortion energy in Fig. ^(c) for NbaSn at i5 = 
(^ = 1) has a small negative value of ^ —4 meV/cell, 
indicating an unstable zone center phonon, since this 
phonon has the same space group as the tetragonally 
distorted A15 cell, but with c = a. The calculation 
for the r2 phonon (pairs of A atoms move toward or 
away from each other with the same displacements for 
all three orthogonal chains) also shows a small, negative 
phonon energy for NbsSn, indicating an instability. Cal- 
culations for the ri2 and r2 phonons for VaSi and VsGe 
also show anomalies. However, for these materials these 
phonons have negative or very small positive energies, 
e.g., their phonon energies are less than 1 meV/cell for 
the ri2 phonon at a displacement of 0.003a, where a is 
the cubic lattice constant. For comparison, zone center 
phonon energies at the same displacement for a typical 
intermetallic compound such as CuaPt in the cubic LI2 
structure are ~ 25 meV/cell (where the unit cell contains 
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FIG. 3. Distortion energies for tetragonally distorted A15 
structure and also see captions of Fig. S. 
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4 atoms). Note that since we estimate the accuracy of 
our LAPW calculations for phonon energies to be approx- 
imately 2 meV/cell, we can argue that the zone center 
phonons and r2 are extremely soft at best, and are 
on the borderline of being stable or unstable. 

First-principles electronic structure calculations have 
shown that the DOS near Ep for V^Si and NbaSn have 
very sharp structures, in qualitative agreement 

with early model calculations of the DOS, and we find 
similar results here. What is particularly interesting, 
however, is that we find explicitly that the DOS singular- 
ities, and the associated Fermi surface nesting features, 
are removed following the sublattice distortions and the 
resultant structural phase transitions. We will discuss 
this more fully in a future publication. 

In summary, it appears that many of the unusual prop- 
erties of the A15 materials are describable in terms of 
"standard" LDA electronic structure calculations when 
sublattice relaxation effects are included. We have shown 
that the A15 materials (VsSi, VaGe, and NbsSn) are in- 
trinsically unstable against shears producing the moduli 
Cii — Ci2 and C44, at T = 0. Our calculations also 
indicate that the zone center phonons and ri2) are 
either unstable or extremely soft at T = 0. Similar stud- 
ies of the full phonon spectrum and the superconducting 
tunneling functions for the A15 materials are being for- 
mulated to verify whether the "standard" LDA methods 
are consistent with all of the experimental data. 

As a word of caution, we urge that elastic stability of 
a "computer-designed" novel material should always be 
checked, as one would not have anticipated the calculated 
instabilities in the A15 materials by examining the total 
energy verses volume curves alone. 

The authors appreciate helpful discussions with M. 
J. Mehl, and gratefully acknowledge the support of the 
Campus Laboratory Collaboration Program of the Uni- 
versity of California. 
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TABLE I. LAPW calculated and experimental elastic constants (in GPa) (from Reference |1|). Calculations were performed 
at the experimental volume using orthorhombic (for Cn — C12) and monoclinic (for C44) distorted cells. "Unrelaxed" refers to 
calculations in which atomic positions are held fixed at their cubic, ideal sites, while "relaxed" refers to calculations in which 
sublattices are allowed to relax. The elastic constants were extracted from a least square fit of distortion energies to the form: 
IS.E — fc2(5^ + kiS"^. "Transforming" ("Non Transforming") denotes an experimental sample that undergoes (does not undergo) 
a cubic to tetragonal phase transition. The "unrelaxed" Cn — C12 values can also be extracted using a tetragonal distortion: 
such extracted values are within 6% of those reported here. 
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